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We study charged brane-world black holes in the model of Randall and Sundrum in which our 
universe is viewed as a domain wall in asymptotically anti-de Sitter space. Such black holes can 
carry two types of "charge" , one arising from the bulk Weyl tensor and one from a gauge field trapped 
on the wall. We use a combination of analytical and numerical techniques to study how these black 
holes behave in the bulk. It has been shown that a Reissner-Nordstrom geometry is induced on the 
wall when only Weyl charge is present. However, we show that such solutions exhibit pathological 
features in the bulk. For more general charged black holes, our results suggest that the extent of 
the horizon in the fifth dimension is usually less than for an uncharged black hole that has the same 
mass or the same horizon radius on the wall. 

PACS numbers: 04.50.+h;98.80.Cq;12.10.-g;11.25.Mj 



I. INTRODUCTION 

In many of the brane-world scenarios, the matter fields 
which we observe are trapped on the brane (see also 
for older proposals). If matter trapped on a brane 
undergoes gravitational collapse then a black hole will 
form. Such a black hole will have a horizon that extends 
into the dimensions transverse to the brane: it will be a 
higher dimensional object. 

Within the context of the second Randall-Sundrum 
(RS) scenario it is important that the induced metric 
on the domain wallQ is, to a good approximation, the so- 
lution predicted by standard General Relativity in four 
dimensions. Otherwise the usual astrophysical properties 
of black holes and stars would not be recovered. 

In a recent paper the gravitational collapse of un- 
charged, non-rotating matter in the second model of RS 



*In this paper, we use the terms 
interchangeably. 



'domain wall" and "brane" 



was investigated. There it was proposed that what would 
appear to be a four-dimensional black hole from the point 
of view of an observer in the brane-world, is really a five- 
dimensional "black cigar" , which extends into the extra 
fifth dimension. If this cigar extends all the way down 
to the anti-de Sitter (AdS) horizon, then we recover the 
metric for a black string in AdS. However, such a black 
string is unstable near the AdS horizon. This instability, 
known as the "Gregory-Laflamme" instability |^] , implies 
that the string will fragment in the region near the AdS 
horizon. However, the solution is stable far from the AdS 
horizon. Thus, one may conclude that the late time solu- 
tion describes an object that looks like the black string far 
from the AdS horizon (so the metric on the domain wall 
is close to Schwarzschild) but has a horizon that closes off 
before reaching the AdS horizon. A similar effect occurs 
when there is more than one extra dimension transverse 
to the brane ||^. These conclusions are supported by an 
exact calculation for a three dimensional RS model |^ . 

In this paper, we consider black holes charged under 
gauge fields which are trapped on the brane. The flux lines 
of such gauge fields can pierce the horizon only where it 
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actually intersects the brane. The bulk theory is the same 
as for the uncharged case so one might expect that the 
black cigar solution would still describe the bulk metric 
of such a charged brane- world black hole. The effect of 
the charge might simply be to modify the position of the 
brane in the bulk spacetime. If this were the case, then 
we might be able to repeat the analysis of ^ by start- 
ing with the black string metric and solving the Israel 
equations appropriate for the presence of a gauge field 
on the brane. However, in the Appendix we prove that 
this is not possible. It is still conceivable that the bulk 
metric is the same as that of the black cigar, but unfor- 
tunately the form of the cigar metric is not known. We 
are therefore forced to study charged brane-world black 
holes numerically. 

A recent paper has claimed to give a solution de- 
scribing a non-charged black hole in the RS scenario. By 
using the brane-world Einstein equations derived in [ pT| , 
it was shown that a Reissner-Nordstrom (RN) geometry 
could arise on the domain wall provided that the bulk 
Weyl tensor take a particular form at the wall. We re- 
gard this solution as unsatisfactory for two reasons. First, 
there is no Maxwell field on the domain wall so the black 
hole cannot be regarded as charged]^. Secondly, only the 
induced metric on the domain wall was given - the bulk 
metric was not discussed. The solution is simply a solu- 
tion to the Hamiltonian constraint of general relativity 
and gives appropriate initial data for evolution into the 
bulk. Until this evolution is performed and boundary 
conditions in the bulk are imposed, it is not clear what 
this solution represents. For example, it might give rise 
to some pathology such as a naked curvature singularity. 
We would then not regard it as a brane-world black hole, 
which should have a regular horizon |^,^ . One aim of the 
present paper is to evolve the initial data of in order 
to understand what this "solution" really describes. 

The second aim of this paper is to study brane-world 
black holes that are charged with respect to a Maxwell 
field on the brane. We start by solving the Hamiltonian 
constraint on the brane to give an induced metric that 
is close to, but not exactly, Reissner-Nordstrom. The 
"charge" of jl^ arises as an integration constant in the 
metric. We then evolve this "initial" data away from the 
domain wall in order to study the resulting bulk space- 
time. Our solution to the Hamiltonian constraint is based 
on a metric ansatz that is almost certainly not obeyed by 
the true solution describing a charged brane-world black 
hole. However, we expect our ansatz to be sufficiently 
close to the true solution that useful results can be ob- 
tained without a knowledge of the exact metric, just as 




t In the AdS/CFT interpretation |T|] of the RS model, this 
black hole must be charged with respect to a (7(1) subgroup 
of the dual CFT. 



Our results suggest that it is more natural to take the 
"charge squared" parameter of to be negative than 
positive since the latter gives an apparent horizon that 
grows relative to the black string as one moves away from 
the brane. For black holes charged with respect to a 
Maxwell field, we find that the horizon shrinks in the fifth 
dimension. In both cases (and for black holes carrying 
both charges), we obtain a numerical upper bound on 
the length of the horizon in the fifth dimension. We find 
that increasing either type of charge tends to decrease 
this length, even if the horizon radius on the brane is 
held fixed. 

It is worth emphasizing that this paper is quite distinct 
from recent papers which have appeared on the subject 
of charged black holes in brane-world scenarios [^-p^. 
This is because these papers all study the effects of bulk 
charges on the brane-world geometry, whereas our anal- 
ysis deals with gauge degrees of freedom that are truly 
localized on the brane. One consistent interpretation of 
the RN solution of would be as the induced metric 
on the brane in the (bulk) charged black string solution 
of ||l^,|l^ • However, in this paper we will study whether 
sense can be made of this solution without introducing 
bulk gauge fields. 

Related numerical work on uncharged brane-world 
black holes has recently appeared in [Q. The differ- 
ence between that paper and the present work is that we 
will prescribe "initial" data on the brane and evolve it in 
the spacelike direction transverse to the brane, whereas 
in [[l6| , initial data was prescribed on a spacelike hyper- 
surface and evolved in a timelike direction. 

The outline of this paper is as follows. First, we set up 
the basic notation and formalism for a covariant treat- 
ment of the second brane-world scenario of Randall and 
Sundrum. Next we solve the Hamiltonian constraint for 
"initial" data on the brane and obtain a RN solution with 
small corrections. We then numerically evolve the solu- 
tion into the bulk subject to the constraint that the met- 
ric solve the vacuum Einstein equations with a negative 
cosmological constant. Finally we discuss the properties 
of the resulting bulk spacetime. 

II. FORMULATION AND STRATEGY 

A. Covariant formulation of brane-world gravity 

We shall be discussing a thin domain wall in a five 
dimensional bulk spacetime. We shall assume that the 
spacetime is symmetric under reflections in the wall. The 
5-dimensional Einstein equation is 

<^''K^-l^'W^R^4^'%., (2.1) 

where = SttGs and G5 is the five dimensional Newton 
constant. The energy-momentum tensor has the form 

= -As^^V- + 5{x)hXh^, + T^,]. (2.2) 
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In the above, the brane is assumed to be located at x = 
where x is a Gaussian normal coordinate to the domain 
wall. X = is the fixed point of the Zi reflection sym- 
metry. A and A denote the bulk cosmological constant 
and the domain wall tension respectively, /i^i, is the in- 



duced metric on the wall, given by /i^^ 



(5) 



where is the unit normal to the wall. The effect of the 
singular source at x = is described by Israel's junction 
condition |17| 



T 



(2.3) 



Here, K^^ denotes the extrinsic curvature of the domain 



wall, defined by if^i, = h^^h'^V pUa- In equation (|2.3|), we 
are calculating the extrinsic curvature on the side of the 
domain wall that the normal point into. This is because 
we want to evolve initial date prescribed on the wall in the 
direction of this normal. Using the Gauss equation and 
the junction condition, we recover the Einstein equation 
on the brane 1 11 1 : 



('^^G^, = -A4 V + SnGiTp, + 4Trp, - E^,, (2.4) 



where 



A4 = 



G4 
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B. Strategy 

We adopt the following procedure: We take a certain 
charged black hole geometry for the brane. When we 
solve for the bulk, we Wick rotate twice. This gives a 
Kaluza-Klein bubble spacetimc ||2l|,^ from which we ob- 
tain boundary conditions at the condition on the bubble 
surface. Wick rotating back gives boundary conditions 
at the bulk horizon for our problem. The Kaluza-Klein 
bubble spacetime is reviewed in Appendix B. 



C. Metric and field equations 

We assume that the induced metric on the brane takes 
the form 



ds^ = ~U{r)dt'^ + 



U{t) 



+ r'^dnl, 



(2.13) 



where dO^I ~ ^9"^ + sin^ ddi^' 



A5 



rri rria 



J.2 



and E^i, is the 'electric' part of the 5-dimensional Weyl 
tensor: 



(2.8) 



^. Note that this is a guess. 
It is unlikely that the exact metric describing a brane- 
world black hole would have precisely this form - in gen- 
eral one would expect the coefficients of dt^ and dr^ to 
be independent (when the coefficient of dn| is fixed as 
, r^). However, we know that the induced metric describ- 
' ing a charged black hole should be close to Rcissncr- 
Nordstrom, which can be written in this form, so our 
■6) ansatz is probably quite a good guess. We expect that 
deviations from the exact metric will give rise to patholo- 
■ 7) gies when this initial data is evolved into the bulk. Even 
so, the analysis of shows that it is possible to extract 
quite a lot of information from a pathological solution. 
The function U{r) will be determined from the Hamil- 
tonian constraint equation below. The bulk metric is 
assumed to take the form 



We shall now specialize to the RS model. This has 
6 



A, 



which implies 



A4 = 0, G4 



G5 



(2.9) 



(2.10) 



The matter on the domain wall will be assumed to be a 
Maxwell field. This implies T = 0, so we can rewrite the 
Einstein equation as 



ds2 ^ N{x,rfdx^ - e2°(x.-)f/(^)rfi2 

„26(x,r)j 2 



Uir) 



(2.14) 



A'' is the lapse function which describes the embedding 
geometry of the hypersurface spanned by (t, r, 9, ip) dur- 
ing the evolution in the bulk spacetime. 

The extrinsic curvature of a hypersurface of constant 
X (with unit normal n = Ndx) is given by 



Kl = ^, K^^ and K'.^K^^^, (2.15) 



b_ 

N 



= 8^G4r^. - ^T^pT^ - E^,. (2.11) 



where a dot denotes d^. The spacetime is described by 
the evolution equation. 



The Israel equation gives the extrinsic curvature of the 
wall: 



(2.12) 



= n(^^^^R^ - KKI^ + D^'D.N, (2.16) 



the Hamiltonian constraint equation. 



(2.17) 
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and the momentum constraint equation, 



D^IO^ - D,K = 0. 



(2.18) 



Here and (^^i? are the Ricci tensor and Ricci scalar 

on hypersurfaces of constant x- 

III. BRANE AND BULK GEOMETRY 

A. Brane Geometry : Charged black hole "initial 
data" 

The action for the MaxweU field on the brane is taken 
to be 



S = 



1 



giving energy-momentum tensor 



T — 



1 



F F P — -h F FP'^ 



(3.1) 



(3.2) 



The field strength F is related to a potential ^4 by F = 
dA. The equations of motion are satisfied if we take 
A = -<^{r)dt with <i>(r) = Qjr. This gives 



T — 



1 



87rG'4 



-diag(U,-U"\r2,r2 sin^ 



(3.3) 



This can be s ubstit uted into the right hand side of the Is- 
rael equation (2.12) to give an expression for the extrinsic 
curvature. This can then be substituted into the Hamil- 



tonian equation (2.17), along with our metric ansatz to 
give an equation for U(r). Solving this equation givefl 



\J(r) = 1 - 



264^/ /? + Q2 ;2g4 



20r6' 



(3.4) 



where M and /3 are arbitrary constants of integration. 
Substituting into the Einstein equation on the domain 
wall gives 



* It is interesting to compare this form for U{r) with the 
behaviour expected from the linear perturbation analysis of 
the second RS model |^,p^|-|2o|. In linearized theory, U(r) = 
1 — (^(r) where 0(r) is the Newtonian potential. For r (., 
the leading order corrections to (/)(r-) are expected to be pro- 
portional to GiMl^ /r'^ and ^^Q^/r**. Such terms are not 
present in our expression for U{r). However, we shall be in- 
terested in black holes for which G4M 2> i, so these correction 
terms will be small compared with terms like {GuM/r)'^ and 
{GaMQ /r^Y , which would be neglected in linearized theory. 
In other words, the RS correction are dominated by post- 
Newtonian corrections so it is not appropriate to compare 
U{r) with the linearized results beyond leading order. 



E„ 



f3 I 



2n4 



2r» 



diag(U,-U~\r2,r2 sin^ ■ 



(3.5) 



It is interesting to compare —E^^i, with SnG^Tfj^iy since 
these quantities appear on an equal footing in the effec- 
tive Einstein equation (2.4). It is clear that the constant 
of integration /? is in some sense analogous to , which 
is why the authors of ||l^ obtained a RN solution. How- 
ever, since their solution did not have a Maxwell field, it 
cannot really be regarded as a charged black hole in the 
usual sense. Rather it carries "tidal" charge associated 
with the bulk Weyl tensor. (3 might be regarded as a five 
dimensional mass parameter. 

We shall only consider initial data that corresponds to 
an object with an event horizon (in the four dimensional 
sense) on the domain wall. In some cases there may be 
more than one horizon. We shall use r_|_ to denote the po- 
sition of the outermost horizon, i.e., the largest solution 
of U{r) = 0. This has to be found numerically except 
when Q — 0. 

Our "initial data" is given by 



a 
N 

c 

N 



b 

N 



i9l 

2 r4 ■■ 



i9l 

2 r4 ■ 



a{x = 0, r) = &(0, r) = c(0, r) = 0. 



(3.6) 

(3.7) 
(3.8) 



We shall study the following cases: 
(i) No electromagnetic charge, i.e., Q = Q. In this 
case, the induced metric on the domain wall is exactly 
RN lliol. The horizon radius is 



M + v/M2 - /3. 



(3.9) 



The induced metric has a regular horizon if /? < M^. 
Note that there is nothing to stop us choosing /? to be 
negative, which emphasizes the difference between the 
solution of |l^ and a charged black hole. If we take (3 to 
be negative then the induced metric has only one horizon, 
instead of the two horizons of a non-extreme RN black 
hole. 

(ii) No tidal charge, i.e., (3 = 0. In this case, the 
induced metric on the domain wall is Reissner-Nordstrom 
with a correction term. Note that —Ep,i, is non-zero but is 
of order 1/r^, which suggests that the total "tidal energy" 
on the wall is zero. 

We shall also consider the general case (ill) Both 
charges non-zero, i.e., /3 7^ 0, Q 7^ 0. 



B. Bulk Geometry 



Th e bu lk geometry is obtained by integrating ( 2.15 ) 
and (2.16) in the ^-direction numerically. We use the 
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standard 'free-evolution' method, that is we do not solve 



the constraint equations ( 2.17 ) and ( 2.18 ) during the evo- 
lution, but instead use them to monitor the accuracy of 
the simulation. 

We obtain the solution numerically in the region r+ < 
r < with Te ^ 5r+ . Boundary conditions at r = r+ 
are specified by first Wick rotating x — ^T, t = ir, 
which takes the metric to a Kaluza-Klein bubble met- 
ric (see Appendix B). Therefore we can apply the nu- 
merical techniques that are used in the study of Kaluza- 
Klein bubbles , although the physics of Kaluza-Klein 
bubbles is unrelated to the physics of black holes. It 
was shown in the Appendix of that at the inner 
boundary r = r+, a and b evolve synchronously, that 
is, a(T, r+) — b{T,r+). Analytically continuing back 
to our original spacetime yields the boundary condition 
a(x, r+) = 6(x, r+). The evolution equation for the trace 
of Kfj^i, and the momentum constraint are also used at 



At the outer boundary r 



we assume the 



components of the extrinsic curvature [equation (2.15)] 
faU off like -!/£ + O (r-*) [cf. (ja^)]. We apply the 
geodesic gauge condition (slicing condition), iV = 1. 

We use the Crank-Nicholson integrating scheme with 
two iterations [ p5| . The numerical code passed conver- 
gence tests, and the results shown in this paper are all 
obtained to acceptable accuracy. 

We were only able to solve numerically in a region near 
the domain wall with a maximum value for x oi 0{1). 
This is because the volume element of surfaces of constant 
X decreases exponentially as one moves away from the 
wall, just as in pure AdS. The evolution was stopped 
when y/—g became too small to monitor accurately. 

We are interested in how charge affects the shape of the 
horizon, in particular how far it extends into the fifth 
dimension. This will be measured by the ratio of the 
physical size of the apparent horizon r+e'^''-^'''+-', to that 
of a black string with the same horizon radius r+ on 
the wall^. The size of the black string apparent horizon 
in the bulk is r+e"-^/^, so the ratio is 



Rix) 



„c(x,r+)+x/ 



(3.10) 



We remark that the only apparent horizon that appears 
during the x^^volution is at r = r_|_. Here we define 
apparent horizon as the outermost region of negative ex- 
pansion of the outgoing null geodesic congruences, where 
we define the expansion rate, as 



_(3) 



Kab, 



(3.11) 



^ The reason for measuring the size of the horizon relative 
to that of the black string is because we want to distinguish 
the closing-off of the horizon from the exponential collapse of 
hypersurfaces of constant x arising from the AdS nature of 
the geometry. 



where s" = {1/ ^grr)dr is an outwards pointin g uni t vec- 
tor in the 3-dimensional metric. We checked ( ^.11 ) dur- 
ing the evolution and confirmed its positivity for r > ?■+ . 

Our initial conditions give the behaviour of the ratio 
R{x) near the brane: 



^(X)lx=o = -^^<0, 



and 



^(x)lx=o = 



3g2+/3 £2g4 



2rl 



(3.12) 



(3.13) 



For model (i) (Q = 0), i?(x)lx=o - 0, but R{x)\x=o = 
l3/r\. This gives R{x)\x=o < for the case with 
/3 < 0, which indicates that the ratio decreases, while 
i?(x)lx=o > for the case with /3 > 0, which indicates 
that the ratio increases. We have plotted the numerical 
results for this ratio in Fig.^ (a) and (b) (henceforth we 
shall set £ = G4 = 1 and assume M 1, as appropri- 
ate for an astrophysical black hole.). Fig.|l| (a) and (b) 
suggests that a negative value for (3 is the natural choice 
since the apparent horizon grows (relative to the black 
string) in the fifth dimension when (3 is positive. 

For model (ii) (/? = 0), Rix)\x=o < and the ra- 
tio always decreases [see Fig.^ (c)]. Model (iii) (Q 7^ 
and /3 7^ 0) is non-trivial. We present numerical re- 
sults in Fig.|. The plot is for M = 5, Q = 3 and 
(3 — 0, ±5, ±10, ±15, where (3 — \b is close to the ex- 
tremep^ case for this choice of Q. The qualitative fea- 
tures are combinations of the plots in Fig.|l|. Note that (3 
seems to have the greatest effect on the bulk evolution. 
Again, the case with negative /? appears to be the natural 
choice since positive (3 gives a growing horizon. 



C. Bulk Geometry: extent of the horizon 

In this section, we shall estimate how far the horizon 
extends into the fifth dimension by combining analytical 
and numerical work. Following the conjugate points the- 
orem 26 1, we shall show that for a charged black hole. 



the trace of the extrinsic curvature diverges at a finite 
distance from the brane. The trace of the evolutional 
equation is given by 

k^^^^R-K^ + ^ = -K^,K^^'' + (3.14) 

where we used the Hamiltonian constraint in the second 
line. Now define k^i, as 



** By extreme, we mean that U{r) has a double root at 
r = r+. 
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— : — + fc^. 



(3.15) 



The trace part of fc^,y, k — k^, is expected to measure 
the volume expansion relative to the AdS "background" 
geometry. In term of fc^,y, (3.14) can be written as 



k~-k + -k' 
£ 4 



kiiuk^ ^ 0, 



(3.16) 



where fc^^, is the traceless part of fc^^. On the brane the 
"initial" condition is 



which implies 

For the case with Q 0, 

k^i/k^ > 0, 

so 

fe|branc 

This implies that there is a x = xo such that 
fc = fco < 



From ( 3.16 ), one obtains 



m 



< 1 



^Ifcol 



o2(xo-x)/^ 



(3.17) 
(3.18) 
(3.19) 
(3.20) 
(3.21) 

(3.22) 



from which it follows that k diverges before x — Xcrit, 
where 



^1 / 

Xcr2t =X0 + 



e\ko\ 



(3.23) 



The divergence in k implies that K also diverges. Near 
X = Xcnt, |fc| behaves like 



k < 



Xcrit — X 



(3.24) 



The case with Q = is more diffic ult to analyze because 
fc|brano = 0. Wc Can usc equation ( 2.15| ) (with iV = 1) to 
give 



k I J- J- _|_ — 2 



2x 



(3.25) 



where we have used the synchronous evolution boundary 
condition a — b a± r ^ r^. In Fig.^, we have plotted 
a + c + 2x/i at r = r+. It is clear from this plot that k 
becomes negative in the bulk when (3 < 0. In fact k also 
becomes negative when /? > 0. Thus, even in the Q — 
case, there exists a x = Xo such that k = kg < 0. The 
above argument can then be used to show that when 



(5 = and (3 ^ Q, K must diverge before x = Xcrit, 
where Xcrit is given by equation ( 3.23 ). We have therefore 
proved that if Q 7^ or /3 then K diverges before 

X ~ Xcrit- 



It follows from equations (3.24) and (B.25) that 



{a + c)\r=r^ < 21og(xc 



(3.26) 



/— g tends to zero at least as fast as 



which implies that 

iXcrit - X)"* as X ^ Xcrit- 

Conservatively, the divergence of K indicates that the 
geodesic slicing has broken down (when N = 1, d-^ is the 
tangent vector of spacelike geodesies), in other words a 
caustic has occurred. The numerical study therefore can- 
not be extended further using this slicing. This has, how- 
ever, a physical meaning because the apparent horizon is 
located at constant r = r-|_ in the bulk. The horizon will 
encounter the caustic before reaching the AdS Cauchy 
horizon. The caustic can therefore be viewed as the end- 
point of the horizon, i.e., the tip of the black cigar. Our 
analysis has only shown that the geodesic slicing must 
break down at the caustic so, in principle, this point may 
be regular, i.e., there may exist a coordinate chart that 
covers a neighbourhood of this point^ However, we do 
not regard this as very likely. Our guess for the induced 
metric on the domain wall is unlikely to be exactly cor- 
rect, so in general we would expect some pathology such 
as a naked curvature singularity to appear in the bulk. 
We cannot check whether curvature invariants diverge at 
X — Xcrit since our numerical evolution cannot be ex- 
tended as far as x = Xcrit- 

Wh ether the bulk solution is regular or not, equation 
(|3.23 ) gives us an upper bound on the extent of the hori- 
zon in the direction transverse to the domain wall, i.e., 
the length of the black cigar. We have plotted this up- 
per bound in Fig.^ taking the values for xo and fco at 
the endpoint of our numerical evolution. The first graph 
shows how Xcrit depends on Q and (3 when M is fixed. 
Note that when Q = — 0, the numerical solution is 
simply the black strin^^, which has Xcrit = 00. Increas- 
ing Q clearly has the effect of decreasing Xcrit , which is 
not surprising since increasing Q also shrinks the horizon 
radius on the domain wall r^. Perhaps more surprising 
is that making (3 more negative also appears to decrease 
Xcrit even though t his increases the horizon radius on the 
wall [see equation (3.E)]. The solid curve on this diagram 
has both AI and r+ fixed. It is clear that Xcrit decreases 
along this curve as Q or /5 increases. 



^^It is not even clear from our analysis whether the caustic 
occurs at a single point or is spread over a region of spacetime. 

^"'"The reader may find this surprising since the black string 
is unstable Q, and small numerical errors might be expected 
to act as perturbations. However, the string is unstable to 
long wavelength perturbations, and the numerical errors will 
only be relevant at short wavelengths. 
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The second graph of Fig.^ plots the same curve (fixed 
M and fixed r+) for different values of M . The trend 
seems to be the same in each case. 

The final graph of Fig.^ is for fixed r+ (rather than 
fixed M). Increasing (3 appears to decrease Xcrit when 
Q is small but has no significant effect when Q is large. 
When /? is non-zero, increasing Q has the effect of ini- 
tially slightly increasing Xcriti but ultimately decreases 
it substantially. The gross trend appears to be that in- 
creasing either type of charge leads to a decrease in the 
length of the horizon. 

In most of these graphs, Xcrit < r^^so the extent of the 
horizon in the fifth dimension is smaller than the horizon 
radius on the domain wall, just as for the uncharged black 
cigar. 

IV. SUMMARY AND DISCUSSION 

In this paper we have studied charged black holes in 
the second RS model. We have seen that two types of 
charge can arise on the brane, one coming from the bulk 
Weyl tensor and one from a Maxwell field trapped 
on the brane. Starting from an ansatz for the induced 
metric on the brane, we have solved the constraint equa- 
tions of 4-1-1 dimensional gravity to find metrics describ- 
ing charged brane-world black holes. In the absence of 
Maxwell charge, one can obtain a Reissner-Nordstrom 
solution . If Maxwell charge is included then one can 
obtain a geometry that is Reissner-Nordstrom with small 
corrections. 

Using these induced metrics as "initial" data, we have 
solved the bulk field equations numerically. We have 
found that the RN solution of [0 has an apparent hori- 
zon that grows (relative to the black string apparent hori- 
zon) in the dimension transverse to the brane unless the 
"charge squared" parameter (3 is taken to be negative^. 
It therefore seems unlikely that this solution really cor- 
responds to a charged brane-world black hole. Of course, 
if a bulk gauge field is included then the work of px| ] 
(with /3 > 0) has a natural interpretation as the induced 
metric on the brane arising from the charged black string 
solution of 

If /? < and/or Q ^ then the horizon shrinks rela- 
tive to the black string horizon. For all cases (including 
(3 > 0), we have found that the trace of the extrinsic 
curvature diverges at a finite distance from the brane, 
with the volume element of the spacetime tending to 



the evolution of Kaluza- Klein bubbles was studied 
numerically and it was found that even though negative mass 
bubbles start off with accelerating expansion j23|, the accel- 
eration ultimately becomes negative. It is conceivable that 
something analogous could happen here but we have found 
no evidence for such behaviour. 



zero. For /3 < 0, we have interpreted this as the end 
point of the horizon of the black hole. Our results sug- 
gest that increasing the charges of a brane-world black 
hole will decrease the length of its horizon in the fifth 
dimension, even when the horizon radius on the brane 
is kept fixed. This implies that, by adjusting Q, one 
can change the five dimensional horizon area while keep- 
ing the four dimensional horizon area fixed. One might 
think that this would lead to a difference between the en- 
tropies calculated from these horizon areas, which would 
be bad news for hopes of recovering General Relativity 
as the effective four dimensional theory of gravity on the 
brane. However, the exponential decrease in the volume 
element as one moves away from the brane implies that 
the dominant contribution to the five dimensional area 
comes from the region of the horizon that is closest to 
the brane |^ . Changes near the other end of the horizon 
give only subleading corrections to the five dimensional 
area, allowing the four and five dimensional entropies to 
agree at leading order. 

We suspect that our solutions will generically have a 
curvature singularity at the point where the trace of the 
extrinsic curvature diverges. This is because it seems 
rather improbable that our ansatz for the induced metric 
on the brane should turn out to be exactly right. How- 
ever, we expect that for each value of Q there will be 
some value of (3 for which a small change in our initial 
data would smooth out this singularity, leading to a regu- 
lar geometry describing a brane-world black hole carrying 
Maxwell charge Q. This smoothing would probably not 
significantly affect the position of the "tip" of the hori- 
zon, for which we have obtained an upper bound on the 
distance from the brane. This is to be contrasted with 
the uncharged case in which one takes the induced met- 
ric on the brane to be Schwarzschild. Evolving this into 
the bulk gives the black string metric, for which the sin- 
gularity occurs at the AdS horizon, which is at infinite 
proper distance from the brane along spacelike geodesies. 
A small perturbation of the metric on the brane takes one 
from the black string to the black cigar, which has a reg- 
ular AdS horizon and a black hole horizon with a tip at 
finite distance from the brane. 

For the black string, the stability analysis of 
shows that the horizon extends a distance of order d = 
t\og{GiM/t) into the fifth dimension, so d ^ r+. Our 
results give only an upper bound for d in the charged 
case. It would be nice if the stability analysis could be 
extended to the charged case. However, the instability 
only sets in when the proper radius of the horizon be- 
comes smaller than the anti-de Sitter length scale and 
we were not able to extend our numerical evolution this 
far. Our upper bound seems rather on the large side, 
since it appears to give d r^ for small Q and p. How- 
ever, for large Q and/or /3, figure^(c) shows that d <C r+, 
so our upper bound is probably tighter in this case. 

The main outstanding problem remains to find the ex- 
act bulk metric that describes a brane-world black hole. 
This was solved for uncharged black holes in the 3 di- 
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mensional RS model by using the 4 dimensional AdS 
C- metric in the bulk p. Unfortunately, the higher di- 
mensional generalization of this metric is not known. It 
would be interesting to see whether charged black holes 
in the 3 dimensional RS model could be constructed by 
using the same bulk as in Q but simply slicing along 
a different hypersurface. It would also be interesting to 
use the methods of [p"8|-p0| to find linearized solutions 
describing spherical distributions of matter charged with 
respect to a brane-world gauge field. 
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APPENDIX A: BRANE BENDING AND THE 
BLACK STRING 



where the electric gauge potential has the form 

A = -^{r)dt (A3) 

so that 



F = <^'{r)dt A dr 



(A4) 



where ' denotes differentiation with respect to r. 

Now, as a first guess we might try to support this 
stress-energy on the brane by allowing the brane to bend 
in the black string background in such a way that the ex- 
trinsic curvatures would still satisfy the Israel equations. 

In other words, we allow the position z of the brane 
to depend on the radial direction r. Solving the Maxwell 
equations yields 



(A5) 



To compute the extrinsic curvature of the timelike hy- 
persurface swept out by z = z{r), we introduce an or- 
thonormal basis which consists of a unit normal vector 



(dz — z' dr) , 



where e = ±1, a unit timelike tangent 
and the spacelike tangents 



t : 




(A6) 



(A7) 



(A8) 



One candidate for a black hole formed by gravitational z 
collapse of charged brane-world matter on a domain wall ~ £r sin 9 'c 

in AdS is the black string solution in AdS, which has the 
metric 

ds^ ^ ^{-Uir)dt^ + U{r)-^dr^ +r^dnl + dz^) (Al) ~ £r 89 



(A9) 



(AlO) 



It follows that the non- vanishing components of the ex- 
where f7(r) = 1 - 2GiM/r. As discussed in g, surfaces trinsic curvature in this basis are 
of constant z trivially satisfy the Israel matching condi- 
tions provided that the tension satisfies A — ±S>Ik^I. — f _|_ If/'^z'), (All) 

Thus, we may slice the spacetime along such a surface £\/l + lIz^ ^ 

of constant z, and match to a mirror image, in order to 
obtain the Schwarzschild solution on the domain wall. 

We now want to consider what happens when we allow j^^^ _ j^^^ _ ~^ _|_ j (A12) 



r 



the black hole to be electrically charged with respect to £\/l + lJz^ 

some U{1) gauge field living on the brane. Thus, we must 
add in an extra term to the brane-world stress energy 

tensor of the form £^ I , ,'2 , rr-i i ^ 

J^tt — 7777 \ zz + z -j- u + 



^--4^(^-^^^-4'^-^-^^") (^2) (AI3) 
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Under the assu mption of Zi symmetry, the Israel equa- 
tions reduce to ( 2.12| ). The three independent compo- 
nents of K^^n give three independent equations: 



Ktt = 



K — 



(A14) 



It is straightforward to show that it is impossible to solve 
these three equations simultaneously unless one takes 
(5 = and z =constant, which is the uncharged solu- 
tion of 1^. It is therefore not possible to support the 
stress-energy of a point charge by simply allowing the 
brane to bend in the black string background. It follows 
that the hulk has to change once the brane-world charge 
is introduced. In other words, brane-world charge will in- 
duce changes in the bulk Weyl tensor, and this is exactly 
what we have found in our numerical analysis. 



APPENDIX B: KALUZA-KLEIN BUBBLE 

The double Wick rotation(x ^ it,t ^ it) of the metric 
of Eq. ( 2.13 ) gives us the Euclidean induced metric: 



ds^ = U{r)dT^ + 



W) 



(Bl) 



The largest r = r+ such that U{r+) = is interpreted as 
the position of the bubble surface. Around r = r-^., the 
metric can be expanded 



ds 



U\r+){r^r+)dT^ + 



dr"^ 



U'ir+){r-r+) 
In term of the new coordinate R := — 



r\dVL\. (B2) 



ds^ 



C/'(r+) 



R'-d 



U'{r+)T\-2 



dR^ 



[dnj. (B3) 



We can see easily that the metric will be regular if 
we assume that the r direction is periodic with period 
47r/f/'(r+). 

In the case of U{r) ~ 1 — r^/r^ with A = A = 0, the 
exact five dimensional solution for time-symmetric initial 
data (if^jy = 0) is 
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dsj = 



This is the Witten- bubble spacetime Another ex- 

ample of initial data for a Kaluza-Klein bubble spacetime 
was given in Ref. |2^ and its classical time evolution has 
been investigated in Ref. P3|,p4|. 
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FIG. 2. Ratio of physical size of apparent hori zon to 
size of black string apparent horizon, -R(x) [cf. eq. (3.1C)], 
for nonzero Q and /3. We have set M = 5, Q = 3 and 
/3 = 0, ±5, ±10, ±15 for this plot. The main features are a 
combination of plots in Fig.^. 
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FIG. 1. Ratio of the physical size of apparent horizon to 
size of black string apparent horizon, R(x) [cf. eq. (3.1C)], 
plotted as a function of x- Fig- (a) is for model (i). Lines 
are of /3 = 0, ±0.2A^^ ±0.5M^ ±0.9A/^ and ±A/^ where 
M = 10.0. We see that the qualitative behaviour of R(x) 
depends on the sign of (3. Fig.(b) is for the extremal case of 
model (i), (3 = with different values of M. Results are 
also plotted for f3 — —M^. Fig.(c) is for model (ii) for which 
R{x) is monotonically decreasing. 
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FIG. 3. The quantity a + c + 2x/£ at r = r+ plotted for 
Af = 5, = and /3 = 0, -5, -10, -25 and -50. 
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FIG. 4. Critical value Xcrit [eq.(3.23)] plotted for non-zero 
Q and /3(< 0) black holes, (a) [above left] We have set Af = 5 
for this plot. Note that for the uncharged case, Xcnt = oo. 
The solid curve is for the special cases with r+ = 10.0. (b) 
[above right] Critical value x^rit for combinations of param- 
eters (Q, /3) which produce a black hole with r+ = 10.0. 
M = 5.0, 6.25, 7.5, 8.75 and 10.0 are chosen for these plots. 
The black dots denote the ends of the lines ai l3 — and the 
other ends projected onto the Xcrit = 3 plane, (c) [below] The 
same plot as for (b), but M is specified so as to fix r+ = 10 
for given 11 



